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Abstract 

Functionally  Graded  Solid  Rocket  Propellants  are  being  developed  at  NSWC-Indian  Head  in 
conjunction  with  the  University  of  Maryland.  The  approach  being  used  treats  these  propellants  as 
typical  Functionally  Graded  Materials  (FGMs),  which  by  definition  are  structures  that  possess 
gradual  variations  in  material  behavior  that  enhance  material  and/or  structural  performance.  For 
functionally  graded  propellants,  Twin  Screw  Extrusion  (TSE)  processing  is  used  to  continuously 
vary  the  composition  throughout  a  grain  in  a  controlled  manner.  As  a  result,  TSE  processing 
allows  the  burning  rates  of  propellants  to  be  tailored  as  a  function  of  burning  web  thickness.  This 
in  turn  will  allow  for  direct  Thrust  Magnitude  Control  (TMC)  for  a  solid  rocket  motor,  which  has 
proven  difficult  to  achieve  in  the  past. 

To  realize  the  benefits  of  functionally  graded  propellants  in  rocket  motors,  an  Inverse  Design 
Procedure  (IDP)  is  being  developed  that  couples  processing,  property,  and  performance  models 
with  mathematical  optimization  techniques  to  enable  designers  to  determine  realistic  gradient 
architectures  that  will  meet  the  performance  objectives.  An  essential  part  of  this  program  is  the 
development  of  a  model  that  is  capable  of  predicting  the  ballistic  performance  of  functionally 
graded  propellants.  The  development  of  such  a  model  requires  that  new  parameters  be  taken 
into  account  that  would  not  be  considered  for  a  conventional  solid  propellant.  For  instance,  not 
only  will  the  burning  rate  change  as  a  function  of  pressure  and  initial  temperature,  but  the  burning 
rate  characteristics  will  also  change  as  a  function  of  position  in  the  grain.  Furthermore,  the 
density  and  the  thermochemistry  associated  with  the  graded  architecture  of  the  propellant  will  be 
changing  in  a  continuous  manner,  whereas  the  grain  will  conventionally  be  represented  by  a  web 
with  discrete  thickness.  The  details  of  the  ballistics  model  and  the  effect  on  performance 
predictions  related  to  discrete  representations  of  continuously  graded  architectures  will  be 
discussed  in  the  paper. 


I.  Introduction 

Solid  propellants  have  provided  one  of  the  main  propulsion  sources  for  rockets  of  all 
types  over  the  past  60  years.  With  all  of  the  improvements  that  have  been  made  in  rocketry  in 
the  past  century,  solid  rockets  have  reached  a  level  of  reliability  to  enable  their  routine  usage1.  In 
recent  years  the  focus  in  the  solid  propellant  industry  has  been  to  increase  the  performance  of 
the  propellants,  improve  their  insensitive  munitions  (IM)  characteristics,  as  well  as  to  address 
some  environmental  concerns  that  come  about  from  the  toxicity  of  many  rocket  propellants’ 
combustion  products.  While  solid  propellant  rockets  have  proven  to  be  relatively  simple  and 
reliable,  one  of  their  major  drawbacks  is  the  inability  to  throttle.  Current  solid  propellants  are 
limited  to  their  own  set  of  burning  rate  characteristics,  which  when  combined  with  a  certain  grain 
geometry  and  nozzle  configuration  limits  the  possible  thrust  profiles  achievable  for  a  given  rocket 
motor.  Some  work  has  been  done  with  pintle  nozzles  to  achieve  some  throttling  capability, 
however  this  adds  complexity  and  weight  to  what  is  ordinarily  a  relatively  simple  design2,3. 
Functionally  Graded  Solid  Propellants  (FGSPs)  have  the  potential  to  change  the  entire  landscape 
of  solid  rocket  motor  design. 
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By  their  very  nature,  FGSPs  can  provide  an  ability  to  create  a  whole  new  set  of  potential 
thrust  profiles  given  very  simple  grain  geometries.  FGSPs  can  offer  many  advantages  never 
before  offered  to  the  solid  rocket  motor  community.  Now  very  simple  grain  geometries  can 
theoretically  be  used  to  obtain  any  thrust  profile  imaginable.  For  instance,  an  end-burning 
propellant  grain  configuration  can  be  created  to  have  a  neutral,  progressive,  regressive,  or  any 
variance  desired.  A  perforated  grain  can  also  be  created  to  have  any  type  of  thrust  profile  as 
well.  By  using  simpler  grain  geometries,  the  stress  concentrations  present  in  more  complex 
rocket  motor  designs  can  be  mitigated,  reducing  the  likelihood  of  structural  failures.  Furthermore, 
FGSPs  can  offer  the  opportunity  to  improve  propellant  mass  fraction  by  allowing  configurations 
such  as  an  end-burner  to  be  used  for  a  mission  that  it  could  not  have  been  used  before.  A  single 
propellant  grain  could  theoretically  act  as  both  booster  and  sustainer  in  systems  that  require  this 
type  of  performance.  This  is  accomplished  by  intentionally  varying  the  propellant  composition  in 
a  controlled  manner  as  a  function  of  web  thickness. 


II.  Objective 

The  primary  objective  of  this  effort  was  to  develop  a  ballistics  model  to  predict  the 
performance  of  rocket  motors  utilizing  FGSPs.  Initially,  the  model  is  designed  as  a  simplistic 
ballistics  model  capable  of  predicting  the  pressure  and  thrust  generated  by  a  rocket  motor  from  a 
given  set  of  input  parameters.  The  primary  input  parameters  are  propellant  composition  and 
burning  rate  as  a  function  of  web  thickness.  The  model  is  designed  to  work  in  unison  with 
property  and  processing  models  for  the  purpose  of  establishing  an  Inverse  Design  Procedure 
(IDP).  The  IDP  couples  processing,  property,  and  performance  models  with  mathematical 
optimization  techniques  to  enable  designers  to  determine  realistic  gradient  architectures  that  will 
meet  their  performance  objectives.  One  of  the  property  models  being  incorporated  into  the 
inverse  design  procedure  is  the  Petite  Ensemble  Model  (PEM)4,  which  is  used  to  attempt  to 
predict  burning  rates  of  the  propellants  in  an  a  priori  fashion.  These  two  models  will  be  coupled 
with  a  Residence  Distribution  (RD)  model  that  predicts  the  gradient  architectures  that  can  be 
processed  using  Twin  Screw  Extrusion  (TSE)  processing5.  Once  completed,  the  IDP  allows  the 
designer  to  start  from  the  performance  needs  of  a  mission,  and  allow  the  IDP  to  calculate  the 
required  gradient  architecture  of  the  propellant  grain  to  meet  the  performance  needs  of  the 
mission  as  well  as  the  processing  conditions  required  for  its  fabrication.  Therefore,  it  was  critical 
to  develop  a  ballistic  model  that  can  simply  and  accurately  predict  the  performance  of  gradient 
architectures  that  will  be  encountered  in  the  IDP. 


III.  Computational  Model 

III. a  Modeling  of  Bum  Rate  Properties 

For  determining  the  material  gradient  and  distribution  in  a  propellant,  a  model  is  required  to 
describe  the  combustion  process.  Based  on  the  prediction  of  this  model,  the  performance  of  the 
FGCEM  can  be  optimized.  For  this  investigation  a  steady-state  combustion  model  called  the 
Petite  Ensemble  Model  (PEM)  is  used4.  The  PEM  is  based  on  a  statistical  treatment  of  the 
propellant  surface  with  multiple  flame  structure  centered  about  characteristic  oxidizer  particles. 
This  model  can  be  summarized  by  the  following  equations: 
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where  Fd  is  the  overall  oxidizer  distribution  function,  D0  is  the  oxidizer  particle  diamter,  Do  is  the 

mean  oxidizer  particle  diameter,  cr  is  the  oxidizer  distribution  function  mode  width  parameter,  r  is 
the  composite  propellant  mean  burning  rate,  rd  is  the  burn  rate  for  a  pseudopropellant,  ad  is  the 
pseudopropellant  oxidizer  mass  fraction,  Rp  is  the  composite  propellant  pressure  coupled 
response  function,  Rpd  is  the  pseduopropellant  pressure  coupled  response  function,  Rv  is  the 
composite  propellant  velocity  coupled  response  function,  and  Rvd  is  the  pseudopropellant  velocity 
coupled  response  function.  PEM  predictions  of  the  variation  in  burning  rate  with  composition  for 
TSE  processed  energetic  materials  can  be  seen  in  Figure  1. 


1.8', 


APfine  mass  fraction  0  2  012  binder  mass  fraction 


Figure  1.  PEM  prediction  of  variation  in  burning  rate  with  composition  for  TSE  processed  energetic  materials 

Mb.  Modeling  of  Ballistic  Performance 

There  are  several  challenges  when  one  considers  the  ballistics  of  rocket  motor  using  functionally 
graded  propellants.  First,  the  burning  rate  parameters  are  continuously  changing.  Preliminary 
results  from  the  previous  section  indicate  that  burning  rates  may  be  intentionally  altered  by  as 
much  as  50%.  Second,  the  thermochemistry  of  the  propellant  is  also  varying  constantly  as  the 
composition  changes.  Finally  the  solid  propellant  density  changes  as  the  gradient  architecture 
changes  throughout  the  web  of  the  propellant.  Codes  currently  available  to  industry  cannot  take 
into  account  all  of  these  challenges. 

As  a  first  attempt  to  predict  the  performance  of  solid  rocket  motors  using  functionally  graded 
propellants,  a  computer  code  as  been  written  that  tries  to  take  into  account  all  of  the  problems 
that  codes  that  predict  performance  of  conventional  solid  propellants  do  not.  Since  this  ballistics 
code  is  to  be  used  as  part  of  larger  IDP,  care  must  be  taken  to  balance  the  accuracy  of  the 
results  with  the  overall  run-time  of  the  code.  Since  the  IDP  may  require  many  iterations,  the  run 
time  of  the  ballistics  code  must  be  kept  to  a  minimum  while  maintain  reasonable  accuracy. 
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Once  the  propellant  gradients  have  been  specified,  the  code  predicts  the  performance  of  the 
rocket  motor.  The  model  assumes  that  each  web  gradient  has  unique  burning  rate  properties, 
thermochemistry,  and  density.  The  code  calculates  the  density  of  the  propellant  as  a  function  of 
web  thickness,  runs  the  NASA  CEA  ’00  computer  code6  to  determine  the  thermochemistry  at 
every  single  web  gradient  that  was  specified,  and  extracts  all  of  the  data  that  is  of  interest  for 
ballistic  calculations,  such  as  the  flame  temperature,  c*,  molecular  weight  of  the  gases,  specific 
heat  ratio,  and  gas  constant.  The  code  creates  tables  of  these  parameters  versus  burning  web 
thickness. 


Once  the  solid  propellant  density,  and  thermochemical  calculations  have  been  made  the  code  is 
ready  to  handle  the  internal  ballistics  and  performance  calculations.  During  the  ballistic 
calculations,  the  code  references  the  tables  created  earlier  as  it  monitors  the  propellant  web 
burned.  It  marches  in  time,  calculates  web  burned,  and  compares  the  web  burned  to  the 
specified  gradients.  If  the  web  burned  is  somewhere  in  between  the  gradients  the  code  will 
linearly  interpolate  for  the  thermochemical  properties  of  the  propellant  and  burning  rate 
parameters. 

The  ballistics  routine  predicting  the  performance  of  the  rocket  motors  can  be  calculated  from  first 
principles.  The  control  volume  considered  is  depicted  in  Figure  2.  The  following  assumptions 
were  made  when  constructing  the  model. 


1 .  Heat  transfer  was  negligible. 

2.  The  mass  added  by  the  burning  propellant  is  at  the  adiabatic  flame  temperature  of 
the  propellant  at  the  specified  gradient. 

3.  The  combustion  gases  behave  as  a  perfect  gas. 

4.  Frictional  forces  between  the  propellant  surface  and  combustion  gas  was  negligible. 

5.  Erosive  burning  was  not  considered. 


dm 


Fig.  2  Differential  Control  Volume 
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When  considering  these  assumptions  and  logarithmically  differentiating  the  principle  equations7'8, 
the  resulting  equations  are  obtained. 
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Energy: 
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Mach  Number: 
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In  addition,  the  static  and  stagnation  pressure  are  related  through  isentropic  flow  relations,  which 
when  logarithmically  differentiated  yields: 
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d(M 2) 
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The  rocket  motor  geometry  is  broken  into  axial  elements  of  equal  size.  Equations  (5)  -  (12)  are 
applied  to  each  element  for  determining  the  incremental  property  contribution  from  each  element. 
The  incremental  contribution  is  then  combined  with  the  incoming  properties  as  is  depicted  in 
Figure  2.  This  combination  then  becomes  the  incoming  set  of  properties  for  the  next  element. 

The  incremental  mass  generation  rate  of  combustion  products,  m„  is  related  to  the  density  of  the 
products,  pp,  the  burn  area,  Ab,  and  the  burning  rate,  rb,  as  follows: 

lilg=PPAbrb  01) 

The  dependence  of  the  burning  rate  on  the  chamber  pressure,  Pc,  is  specified  in  the  traditional  St. 
Robert’s  Burning  Rate  Law: 


rb=aPcn  (12) 

where  n  is  the  burning  rate  exponent  and  the  burn  rate  coefficient,  a. 

Once  the  calculation  has  been  performed  for  each  element,  the  final  ballistic  calculations  can  be 
made,  starting  with  equation  (13).  The  conservation  of  mass  can  be  used  to  determine  the 
change  of  mass,  m,  in  the  combustion  chamber  from  the  exiting  mass  flow  rate,  me ,  and  the  total 

mass  generation  rate  of  combustion  products,  mgt ,  as  follows: 


dm 

-  =  (13) 

The  exiting  mass  flow  rate  is  determined  from  the  exit  throat  area,  At,  using  the  choked  flow 
assumption: 


e 


PA 

c  * 


m 


(14) 


Proceedings  of  the  40th  JANNAF  Combustion  Meeting  (2003) 


The  equation  of  state  for  an  ideal  gas  can  then  be  used  to  relate  the  mass  to  the  state  variables, 
chamber  pressure,  chamber  volume,  Vc ,  and  chamber  temperature,  Tc,  as  follows: 


m  = 


PV 

c  c 

R  T 

S  C 


(15) 


Differentiating  (15)  and  combining  it  with  (13)  yields: 

P,  dV  V  dP,  PV  dT  PV  dRB 


dm  _  d 
dt  dt 


^  P  V  ^ 

c  v  c 

v¥c; 


RgTc  dt  RgTc  dt  RgT;  dt  RgT.  dt 
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The  instantaneous  chamber  volume,  Vc(t),  is  given  by: 


K(t)  =  vu<u  +  \'-^Ldt 

i  p„ 

dV„  ni, 


Differentiating  (17)  yields: 


dt  p. 


(17) 


(18) 


Substituting  (18)  into  (16)  and  rearranging: 


dpc  _  Rj. ; 


dt 


V„ 


(mg  -me)- 


\m„  -  m ,  I — + 


P  ni  p  dT  P  dRa 


Vc  PP  Tc  dt  Rs  dt 
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Where  chamber  pressure  is  solved  incrementally  by: 


dP 

Pc  =PC  h - -At 

Ci  dt 


(20) 


The  time  step,  zlf.  is  represented  by: 


At  =  —  t 
5 


(21) 


Where  t  is  the  propellant  gas  residence  time  in  the  combustion  chamber  and  can  be 
approximated  by: 


t  = 


PV 


(\-n)[RgTin,) 


(22) 


Both  the  time  derivatives  of  Tc  and  Rg  are  calculated  via  a  one  sided  difference  technique  from 
the  thermochemical  data.  For  example: 


dXi  _  Xt  -X. 
dt  At 


(23) 
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Where  X,  can  represent  state  variables  Tc  or  Rg. 

The  thrust  of  the  rocket  motor,  F,  is  calculated  in  the  following  manner: 


Where: 


And: 


F  =  CfPcAt 


(24) 
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Since  s  is  a  known  variable,  (26)  can  be  used  to  solve  iteratively  for  the  nozzle  exit  plane 
pressure,  Pe.  Once  Pe  is  known,  it  is  fed  into  (25)  to  solve  for  the  thrust  coefficient,  Cf,  from  the 
isentropic  expansion  factor  for  the  combustion  products,  y. 

IV.  Discussion  of  Results 

The  ballistics  code  was  written  in  the  C  programming  language.  When  run  on  a  desktop  PC 
using  2  Pentium  III  1.0  GHz  processors,  typical  run  times  are  on  the  order  of  1  second  for  end¬ 
burning  configurations,  and  2  seconds  for  perforated  configurations  depending  primarily  on 
propellant  grain  web  thickness.  As  the  code  was  written,  the  industry  standard  SPP  '97  computer 
code9  was  used  to  verify  the  results  for  conventional  propellants.  Figure  3  shows  that  very  good 
agreement  was  achieved  with  the  SPP  results.  Typically,  the  predicted  pressures  were  within  1% 
once  a  “steady  state”  was  achieved.  Figure  3  displays  the  agreement  for  an  end-burning 
propellant  grain,  but  similar  agreement  has  also  been  achieved  with  a  center  perforated 
propellant  grain.  The  main  discrepancies  came  during  the  chamber-filling  phases.  These 
discrepancies  are  being  investigated.  However,  in  this  particular  example,  which  is  for  a  relatively 
small  rocket  motor,  it  appears  that  the  chamber-filling  event  occurs  too  slowly  in  the  SPP 
prediction.  One  possible  reason  for  this  is  user  error. 


-  7  - 
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Figure  4  is  an  example  of  the  output  from  the  code  when  a  graded  propellant  is  used  in  an  end¬ 
burning  configuration.  It  shows  the  pressure  and  thrust  history  predictions  for  two  separate  end 
burning  propellant  grains,  one  with  a  conventional  composite  propellant,  and  another  with  a 
graded  propellant.  This  example  shows  how  using  just  one  propellant  configuration  could  create 
a  boost-sustain  system  while  utilizing  a  very  simple  grain  geometry  and  potentially  maximizing 
mass  fraction. 
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Fig.  4  Performance  Predictions  for  Graded  and  Conventional  Propellants  in  an  End- 

Burning  Configuration 


Figure  5  shows  the  output  from  the  ballistics  code  for  both  continuously  graded  propellants  and 
conventional  composite  propellants  when  a  center-perforated  grain  configuration  is  used.  In  both 
Figures  4  and  5,  the  propellant  web  starts  out  with  a  high  burning  rate  and  ends  with  a  low 
burning  rate.  These  high  and  low  burning  rates  were  considered  the  limits  for  this  study,  and  as 
such  calculations  were  performed  for  the  high  and  low  burning  rate  propellants  if  they  were 
considered  conventional  propellants.  The  transition  from  high  to  low  burning  rates  is  done  in  a 
continuous  manner.  Curiously,  in  Figure  5  it  is  noted  that  the  graded  propellant  experiences  a 
“hiccup”  in  which  the  pressure  actually  becomes  higher  than  the  case  of  the  high  burning  rate 
propellant.  This  phenomenon  is  a  result  of  the  gradient  selected  for  this  example  and  the 
changing  burning  surface  area  of  the  center-perforated  grain.  The  results  predicted  by  the 
ballistics  code  demonstrate  that  with  sufficient  computational  resources,  it  may  be  possible  to 
generate  nearly  any  thrust  profile  desired  by  utilizing  functionally  graded  solid  propellants.  This 
ballistics  code,  when  coupled  with  the  IDP  discussed  earlier  will  allow  a  designer,  after  making 
some  basic  sizing  calculations,  to  select  a  mission  thrust  profile  and  input  the  profile  into  the  IDP 
and  have  it  return  the  processing  conditions  necessary  to  generate  the  corresponding  gradient 
architecture. 
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Fig.  5  Performance  Predictions  for  Graded  and  Conventional  Propellants  in  a  Center- 

Perforated  Configuration 


V.  Conclusions 

A  ballistics  code  has  been  written  to  model  the  performance  of  rocket  motors  utilizing  functionally 
graded  solid  rocket  propellants.  The  code  has  been  verified  to  give  satisfactory  accuracy  for 
conventional  solid  propellants  by  comparing  the  results  with  the  industry  standard  SPP  ’97 
computer  code  while  maintaining  reasonable  run  times.  This  ballistics  model  when  coupled  with 
processing,  property,  and  performance  models  with  mathematical  optimization  techniques  will 
form  an  Inverse  Design  Procedure  (IDP)  that  will  allow  designers  to  define  a  performance 
requirement  of  a  mission,  make  some  simple  scaling  calculations,  and  allow  the  IDP  to  calculate 
the  required  gradient  architecture  of  the  propellant  grain  to  meet  the  performance  needs  of  the 
mission.  At  the  time  of  publication,  burn  rate  and  small  scale  rocket  motor  testing  was  imminent. 
Appropriate  modifications  will  be  made  if  necessary  to  ensure  that  the  ballistics  routine  will  give 
results  with  a  satisfactory  level  of  accuracy  to  ensure  success  of  the  IDP. 
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Nomenclature 

Burning  rate  pre-exponential  factor 

Burning  surface  area 

Nozzle  Exit  Area 

Nozzle  Throat  Area 

Thrust  Coefficient 

Constant  pressure  specific  heat 

Characteristic  Exhaust  Velocity 


Rate  of  change  of  mass  in  the  rocket  motor 


Rate  of  change  of  pressure  in  the  rocket  motor 


Rate  of  change  of  gas  constant  in  the  rocket  motor 


Rate  of  change  of  temperature  in  the  rocket  motor 


Rate  of  change  of  free  chamber  volume  in  the  rocket  motor 


Thrust 

Mass  exiting  the  rocket  motor  nozzle 

Incremental  mass  generated  by  burning  propellant 

Total  mass  generation  rate  by  burning  propellant 

Combustion  gas  Mach  number  in  the  rocket  motor  bore 
Burning  rate  exponential  factor 
Ambient  pressure 

Combustion  chamber  static  pressure 
Pressure  at  nozzle  exit  plane 
Stagnation  Pressure  in  Rocket  Motor 
Combustion  gas  constant 
Combustion  chamber  temperature 
Combustion  gas  velocity  in  the  rocket  motor  bore 
Combustion  chamber  free  volume 
Initial  combustion  chamber  free  volume 
Nozzle  exit  plane  to  throat  area  ratio 
Combustion  gas  specific  heat  ratio 
Combustion  gas  density 
Solid  propellant  density 
Time  step 

Propellant  gas  chamber  residence  time 
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